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Abstract 

At the Large Hadron Collider, heavy particles may be produced in pairs close to their 
kinematic threshold. If these particles have strong enough attractive interactions they 
may form bound states. Consequently, the bound states may decay through annihilation 
back into the standard model. Such annihilation decays have the potential to provide 
much information about the bound particles, such as their mass, spin, or charges, in a 
manner completely complementary to standard single particle cascade decays. Many of 
the signatures, such as dijet resonances, will be challenging to find, but may be extremely 
helpful in unraveling the nature of the new physics. In the standard model, the only 
novel annihilation decays would be for toponium; these will be hard to see because 
of the relatively large width of the top quark itself. In models with super symmetry, 
marginally visible annihilation decays may occur for example, from bound states of 
gluinos to dijets or tops. If new particles are bound through forces stronger than QCD, 
annihilation decays may even be the discovery mode for new physics. This paper presents 
various theoretical results about bound states and then addresses the practical question 
of whether any of their annihilation decays can be seen at the LHC. 



1 Introduction 



It will undoubtedly be much easier to conclude that the LHC has found new physics than 
to determine what that new physics is. Many of the most well-motivated models have sig- 
natures such as events with a large number of jets and missing energy. While there is little 
background from the standard model for these signatures, the difficultly of sorting out which 
jet corresponds to which parton, working out the underlying topologies, reconstructing the in- 
termediate masses, and isolating the relevant couplings may be an insurmountable challenge. 
Special cases, for example, with multiple leptons, will make things easier, but determining 
couplings and masses for strongly interacting particles may require the precision of a next- 
generation machine. With no such machine on the horizon, it is therefore imperative to 
consider other ways in which parameters of the underlying Lagrangian may be extracted from 
data at the LHC. In this paper, we consider one often neglected possibility: massive, strongly 
interacting particles that are produced in pairs can form bound states which then annihilate 
back into the standard model. Such annihilation decays are subject to an entirely different set 
of systematics than standard cascade decays, and therefore, if they can be seen, provide the 
possibility of generating entirely complementary information about the underlying physics. 

Bound states exist when the potential between the particles is attractive, with a binding 
energy that is larger than the particle's intrinsic width. If the states are too strongly bound, 
such as through new beyond-the-standard model (BSM) forces, we may never see the con- 
stituent particles on their own, that is to say, the bound states will look just like resonances. 
On the other extreme, if the binding energy is much smaller than the width, the equivalent of 
bound state production and annihilation is simply virtual pair production, with little evidence 
of resonance behavior over the continuum. We are interested in the intermediate regime, where 
the particles sometimes bind and annihilate, and other times have single particle decays, since 
this is when the two complementary types of signatures can be available. Conveniently, this 
happens fairly generically if the particles are bound through the strong force 577(3) of the 
standard model and the single particle decays only go through off-shell intermediate states. In 
particular, it happens for many standard BSM scenarios, where the bound states will typically 
amount to a few percent of the cross section of the pair production processes. 

If the products of annihilation decays can be detected and distinguished from the back- 
ground, they will provide a wealth of information. For example, the location of the peak in the 
invariant mass distribution of the annihilation products can lead to mass measurements more 
precise than those from cascade decays with missing energy. The spins of the possible bound 
states are determined by the properties of the constituent particles. The overall production 
and annihilation rates are also sensitive to the color and spins of the particles. 

We begin in Section [2] with a discussion of bound states of the top quark. Although we 
will conclude that toponium is practically invisible at the LHC, this section introduces some 
of the formalism for studying bound states that we will use in other sections. We proceed, in 
Section [3]to look at beyond-the-standard model examples. While the same basic ideas apply to 
any model of new physics with pair-produced colored particles, to be concrete we will discuss 
mostly bound states of color-adjoint fermions, such as gluinos in the minimal supersymmetric 
standard model (MSSM). When a gluino pair is produced close to threshold, it can form 
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a bound state, called gluinonium [JJ El E|. General features of the gluinonium production 
cross section and annihilation rates, including various perturbative corrections, have been 
studied in [U EJ El El El IS]- Applying these results, we will find that annihilation decays 
of gluinonium are relevant for realistic SUSY spectra such as some of the SPS benchmark 
points [10], as well as for more exotic scenarios such as those with a gluino LSP [UJ [12] or 
split supersymmetry [13]. Some possibly observable final states from annihilation decays of 
gluinonium include dijets and tt pairs. In Section HI we study these decay modes through 
Monte Carlo simulation, and compare them to the relevant backgrounds. In Section [5] we 
discuss bound states of particles with other spins, color representations and charges than the 
gluino. These include squark-(anti)squark and squark-gluino bound states in the MSSM. In 
Section [6] we explore what happens when the particles are bound by a new force that is stronger 
than £77(3)qcd- Section [7] summarizes the various results we find, and concludes. 

2 Wide bound states: toponium 

A reasonable place to start is a review of the bound state formalism for heavy colored par- 
ticles in the standard model. In order to study bound states using perturbation theory, the 
mass of the particles which bind should be larger than Aqcd- Bound states of charm quarks 
(charmonium, such as J/ip) and bottom quarks (bottomonium, such as T) have been well- 
measured and are well-understood theoretically. Bound states of the top-quark, toponia, have 
been studied theoretically as well, but have not yet been seen in colliders. In this section, we 
will review what is known about toponia, and consider the prospects for their observation. We 
will see that the tt invariant mass spectrum does show evidence of binding, but the top quark 
width is too large to allow for any detectable amount of annihilation decays at the LHC. In the 
next sections, we will study more general non standard-model scenarios in which annihilation 
decays are more relevant. 

2.1 Bound state formalism 

The appropriate way to study the formation and decay of bound states is with non-relativistic 
QCD (NRQCD). At high energy, relevant to the production of massive particles at colliders, 
QCD is a perturbative field theory. Therefore, at leading order, the dynamics of a pair of 
colored particles of mass m 3> Aqcd can be described by the single-gluon exchange potential 
which has the form 



The color factor C is obtained by evaluating the product of SU (3) generators. Labeling the 
particles 1 and 2, this product is -TfT 2 a = \({Tf) 2 + (T 2 a ) 2 - (Tf + T 2 a ) 2 ), which leads to 



where C\ and C2 are the quadratic Casimirs for the two particles and C(i2) that for the bound 
state. 
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(1) 



C = - (Ci + C 2 - C(i2)) , 



(2) 
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The coupling constant in Eq. ([TJ is defined as 

a s = a^ciQ 1 ) , (3) 

where ao = (\Ca s m)~ l is the Bohr radius of the bound states. This notation is chosen to 
emphasize that the appropriate scale at which a s should be evaluated in the Schrodinger 
equation is associated with typical momentum transfers relevant for the binding, [i ~ 1 , 
as opposed to the much higher energies \i ~ m which are relevant for the hard scattering 
production and annihilation processes. Since a s will eventually be raised to the third power, 
this can lead to a factor of 2 or more in the overall rate. Note that the color factor C depends 
on the color representation of the bound states and therefore there is a weak dependence in a s 
as well. Other quantities, such as the binding energy, will also depend on the representations, 
and to keep the notation concise, we will generally leave this dependence implicit. 

When the binding particle's intrinsic width T is not negligible compared to the binding 
energy, it must be included in the Schrodinger equation. We will see that this is in fact the case 
for the top quark. The spectral density comprising the bound state resonances is contained 
in the Green's function solving 

G(x,£) = ^(x), (4) 

where E = \/I — 2m is the energy from the scattering event available to the system. The 
bound states are produced from incoming partons with momenta k which correspond to much 
shorter wavelengths than the distance over which the bound state wavefunctions have support. 
Therefore, the production process probes this Green's function at effectively a single point and 
all we really need for the collider physics applications is G(0, E). For example, by the optical 
theorem the partonic production cross section including bound state effects is 

An 

a(s) = —ImG(0,E)a (s), (5) 

where (5 = \/l — 4m 2 /s is the velocity of the two particles in the center-of-mass frame and 
<T is the pair production cross section without including any binding effects. We can also use 
G(0, E) as a propagator for the bound states, for example in computing matrix elements for 
the annihilation processes. 

For zero angular momentum, the solution for the Green's function is [T4"] 



-— + V{r)-{E + iT) 
m 



G(0,E) 



m 
Air 



E + iT 



Ca 9 In 



\C\a t 



m 



E + iT 



E 

n=l 



E n 



-(E + iT) -sign(C)V£n 



(6) 



3 



Here E n = E^/n 2 are the energies of the radial excitations of the ground state whose binding 
energy is given by 

E b = ^C 2 a 2 m . (7) 

This solution encodes the resummation of Coulomb gluons to all orders in a s . Further precision 
can be achieved with the inclusion of higher order corrections to the Coulomb potential (see, 
e.g., [HI [16] and references therein). 

We give further details of the bound states formalism with finite width, including expres- 
sions for the production and decay rates in terms of G(0, E), in Appendix^ The annihilation 
decays will really only be relevant for narrow width, where the formulas simplify. The formulas 
in the narrow width approximation are discussed both in Appendix [A] and in Section [3] of the 
main text. 

2.2 Toponium 

For the case of a tt pair, like for the other quarkonia, the color decomposition 

3<g>3 = l©8 (8) 
shows that toponium can be a color singlet or a color octet. Then 

4 C 'a 1 

C 1 = C F = - , C% = Cf — = — - • (9) 

6 lb 

Thus, the singlet has an attractive potential while the octet does not. The top quark's intrinsic 
width is Y t = 1.3 GeV which is of the same order as its binding energy Eb = \C F oi 2 m t = 1.5 
GeV, so we must include the width in the toponium production. 

The cross section for tt production in the color singlet and color octet configurations is 
shown in the first panel of Figure [TJ Since the octet potential is repulsive, it has no resonance 
behavior. The singlet configuration shows evidence for a resonance, but because of the large 
width of the top, the would-be bound state merges with the tt continuum. The rate for 
annihilation decays of the bound state, the color singlet, is small and has to be scaled by 
a factor of 200 to appear on the same plot. This can be seen analytically in the (poor) 
approximation of a narrow toponium, where the rate of the annihilation process (tt) — > gg 

is [IT] 

64 

T (tt)^99 = g?^™* - 4 x 10~ 3 GeV ~ 2 x 10~ 3 • 2T t . (10) 

Most of the production cross section corresponds to processes that end up in the weak decays 
of the tops, but there are also tiny branching ratios into the various annihilation channels. 
The physical cross sections times branching ratios for it at the LHC (14 TeV) are shown in 
the second panel of Figure [TJ Higher-order QCD corrections to the production rate (but not 
the annihilation processes) have been studied recently in [TSJ [TB], but we do not take them 
into account for our estimates. 

In order to find evidence for binding, it is natural to search for the resonance behavior 
in the invariant mass of the decay products. It is also possible to search for the onset of the 
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Figure 1: Near threshold ti production. The left panel shows the partonic cross section, 
separated by color, with the rate for annihilation decays scaled by 200. The right panel shows 
the physical cross section for ti production at the LHC (14 TeV). The rates for annihilation 
decays into dijets and 77 are shown, scaled by factors of 200 and 20000, respectively. 

continuum contribution which looks like a step in the invariant mass distribution. In any case, 
the challenge is to find a subtle feature with a width of order ~ 1 GeV on top of a smooth 
continuum. In terms of overall rates, the best place to look for evidence of toponium might 
be in the reconstructed invariant mass spectrum from the weak decay products. However, 
because the reconstruction of tops is in itself full of uncertainties, and because of the radiation 
broadening of the hadronic part of the top-decay final states, this is not a promising direction 
to finding a subtle 1 GeV wide feature. The dijet annihilation decay signal is similarly hope- 
less. A slightly more promising channel, at least in terms of resolution, might have been the 
annihilation of toponium into photons, which has branching ratio of (8/9)(a/a s ) 2 ~ 0.005 rel- 
ative to the annihilation into dijets [17]. However, for this decay mode da/dM ~ 0.01 fb/GeV 
which leads to too few signal events, especially compared to the standard model background 
which has da/dM ~ 1 fb/GeV. Similar conclusions have been reached in the limit of narrow 
toponium in the past [TSJ 1201 l2Tj . 

We have seen that there are two main difficulties in finding toponium. The first is that the 
rate for annihilation decays is small, of order 0.2% of the production rate. The second is that 
the large width of the top prevents the resonance from being seen on top of the continuum 
background. Instead of showing a resonance shape, the annihilation decays of toponium into 
dijets produce a signal which merges smoothly with the continuum. Both of these problems 
would be ameliorated if the top had a smaller intrinsic width. Figure [2] shows what happens 
to the production and annihilation decay rates as the top quark width T t is decreased. First 
note from the first panel that the resonance becomes much sharper for smaller width. Not 
only does the bound state become well-defined but its first radial excitations become visible 
too. Thus the distinction from the continuum in any decay mode will be sharper. Second, 
note from the second plot the relative importance of the annihilation decays of the resonance 
compared to the continuum. Of course, the top width is fixed, and known. But there are 



5 




Figure 2: Varying the top width. On the left is the ti color singlet partonic production cross 
section and on the right is the rate for annihilation decays. The lowest curves correspond to 
the real top, as in Figure [U and the others, from bottom to top, have the top width decreased 
by factors of 2, 5 and 10. Of note are both the appearance of the bound states as well-defined 
resonances, and the relative fraction of annihilation decays coming from resonance region. 

beyond the standard model scenarios which have phenomenology similar to a narrow top. 
One example is bound states of gluinos, gluinonium, to which we now turn. 

3 Narrow bound states: Gluinonium 

Having found that toponium is practically impossible to see at the LHC, let us consider in 
what situation bound states of colored particles might be detectable. We saw in Figure |2]that 
the strength of the annihilation decay signal is enhanced when the single particle width is 
decreased. If we can neglect the particle's intrinsic width in the bound state formalism, the 
decay of a resonance can be simply expressed in terms of its wavefunction ?/>(x). The relevant 
rates are then 

2 

r a „n~^|^(0)| 2 ~a^m, E h ~a 2 s m. (11) 

mr 

So r ann /i?5 ~ a 2 s a s <C 1. Thus, bound states in this situation will be very narrow compared 
to the binding energy, like in Figure [2j 

As a canonical example, we will now consider bound states of gluinos, gluinonia. We will 
return to other cases in Sections [5] and [6j Unlike the top whose properties we know, the gluino 
has, minimally, two free parameters: its mass m g and its width Tg. In many scenarios the 
gluino decays to a squark and a quark, with a decay rate scaling like 

T~ g ~a s m g . (12) 

This leads to T g ~ Eb and the gluinonium will be broad like the toponium. However, if 
these two-body decays are kinematically forbidden, which happens when all of the squarks 



6 



are heavier than the gluino, the gluino can only decay through off-shell intermediate states 
(see [221123] for detailed studies of the various possible decay modes in this situation). In this 
case, the decay rate is significantly suppressed. For example, if the dominant decay mode is 
to quarks and a chargino or neutralino through an off-shell squark, the rate will scale like 

a em a s mg fm~ g \ A 
9 ~ 16vrsin 2 ^ \mj ' 1 ' 

In principle, this can be as small as one likes if the squarks are heavy, which enhances the 
fraction of annihilation decays. Since the rate goes as the fourth power of mg/mq, the squarks 
only need to be slightly heavier than the gluinos to saturate the rate for annihilation decays. 

Keep in mind, for the annihilation decays to occur, the gluinos have to bind, which only 
happens near threshold. If the gluinos are produced with relativistic speeds, they still undergo 
single-particle decays rather than annihilating. In the case where the squarks are completely 
decoupled like in split supersymmetry [13] . the relativistic gluinos may even form bound states 
with standard model particles, for example, forming R-hadrons (see, e.g., [HJ [T2] [2H |2S] ) . Such 
exotic phenomenology is largely irrelevant to the cross section for annihilation decays in the 
resonance region, which is the place to look for gluinonia. 



3.1 Narrow gluinonia 

Since we will be interested in the limit Tg, r ann <C Ef,, we can use the narrow width approx- 
imation to simplify the bound-state formalism. In the narrow width limit, the production 
and annihilation of the bound states are determined by finding the energy eigenstates of the 
Schrodinger equation from the potential ([T|). The solutions are then wavefunctions anal- 
ogous to the wavefunctions of the hydrogen atom. For S-waves, these are simply functions of 
the radial coordinate: ?/v( x ) = 2^Rn{ r )-, where n is the radial quantum number. In terms 
of the formalism of Section |2J below the kinematic pair-production threshold, the imaginary 
part of the Green's function reduces to 

ImG(O) = ti \M0)\ 2 S(E + E n ) . (14) 

n 

The short-distance part of the bound state production process is determined by the same 
matrix element .Mo that describes the production of the constituent particles in the continuum. 
The matrix element for binding is given by [26] 127] 

■Mbound = yJv'(O) M , (15) 

where M ~ 2m is the mass of the bound state. Note that, although the continuum production 
rate vanishes at threshold \/§ = 2m, this is due to phase space; the matrix element .Mo a t 
threshold is finite. 

For identical particles, a factor of 1/2 needs to be added in ( ITS"]) and the wave function needs 
to satisfy an appropriate (anti) symmetry condition. For an 5"- wave spin-0 color-symmetric (or 



7 



spin-1 color-antisymmetric) bound state of two gluinos (identical fermions), we should take 
the spatial part as 

V^(x) = -L [^( X ) + V(-x)] = v^x) . (16) 

These factors result in an overall factor of 1/2 in the bound state production and annihilation 
rates. Following the convention in the literature, ^(x) in our expressions will still be normalized 
as in the hydrogen atom, while the extra factor of 1/2 will be included explicitly in the 
pref actor. 

The binding energy, Bohr radius, and the wavefunction at x = for the ground state are 
given by 

C 2 a 2 m 2 2 1 C 3 a 3 s m 3 
E b = — , a = — — , ^0 = — 3 = — . 17 

4 Ca s m naQ 8n 

As before, we have made a notational distinction for the strong coupling constant evaluated 
at scales relevant for the binding a s = a s (dQ 1 ) and the strong coupling scale relevant for the 
production, a s = a s (m). To get a feel for the sizes of the quantities in (fT7|) . we can try C = 3 
and m = 300 GeV, giving E h ~ 11 GeV, l/a ~ 56 GeV and \ip{0)\ 2 ~ (40 GeV) 3 . 

Excited states are much less important than the ground state. The radial excitations have 
wavefunctions which at x = are given by ip n (0) = ip(0)/n 3 ^ 2 , so their decay rates and cross 
sections are suppressed by 1/n 3 . All angular-excited states have ^(0) = 0, so their rates 
depend on the derivative of VK X ) with respect to the radial coordinate. These rates are then 
proportional to \ip'(0)\ 2 /m 2 ~ a 5 s m 3 rather than |-?/;(0)| 2 ~ o7^m 3 [27], which suppresses their 
decay rates and production cross sections by a factor of a 2 compared to the ground states. 
For this reason, we will only consider S'-wave gluinonia. However P-wave gluinonia have been 
mentioned in the context of diffract ive processes in [281 EE]. Note also that because their 
annihilation rates are suppressed, they can be corrected at the leading order by subleading 
terms in NRQCD. For example, in quarkonia, P states can have small admixtures of states 
that are made of a qq pair in an S color-octet state + a gluon [30] . 

The color representation of a pair of gluinos can be one of the following: 

8® 8 = 1s©8 s ©8 a ©10 a ©T0a©27s. (18) 

Here S and A subscripts refer to whether the wavefunction is symmetric or anti-symmetric in 
color. Based on ([2]), the color factors are 

1 3 

C\ = Ca = 3 , C 8g = C 8a = -C*a = - , Cio = Cjq = , C 2 7 = —1 , (19) 

so the potential is attractive for 1, 8§ and 8a, zero for 10 and 10, which are not produced 
at leading order [31], and repulsive for 27. Note that the color factors are raised to the third 
power in the expression for |-0(O)| 2 in 017|) . and this can make a large difference in the cross 
section. It enhances the color singlet production over the color octet by a factor of 8. S'-wave 
bound states exist for all the three attractive representations. Since the gluinos are identical 
spin-i Majorana fermions, the color-symmetric states must be pseudoscalars ( J PC = h ) 

1^0(1)) = 4=**, r^o(8s)) c = \J\d abc (20) 



and the color-antisymmetric state must be a vector ( J PC = 1 ) 

| 3 ^(8 A )> c = i=/ abc , (21) 
where 2S+1 Lj describes the rotational quantum numbers. 
3.2 Annihilation decays 

The gluinonia with different spins will have different annihilation modes, so let us consider the 
pseudoscalar and vector bound states separately. Pseudoscalar gluinonia in either the color 
singlet, 1, or adjoint, 8s, representation will annihilate predominantly into two gluons via the 
diagrams A and A in Figure [31 Including also the normalizations from (j!5p , (jlfip and (1201) , 
the resulting rates are [H [2], [7J 

18ira 2 . . , l9 243 , . , 

T(l gg) = = -r<*ffm s , (22) 



and 



T(8s ^99) = ^§ I^8 S (0)| 2 = 2 ^ s alm g . (23) 



The pseudoscalar gluinonia can also decay into quarks via diagrams B and B in Figure [3j 
The annihilation decay to light qq pairs is suppressed by the quark masses by chirality con- 
siderations. However, decays to tt may have a reasonable rate [51 [9]. The branching ratio for 
the 1 gluinonium is [9] 



T(l ->■ tt) 16 / m\ 



1 2 

m t m g 



T(l-+gg) 27 y mf ( m | + m | _ m ^ 
assuming that both stops have the same mass m t ~. For the 8s gluinonium we similarly find 



(24) 



r(8 s -> tt) 20 / m\ m 2 m 



r(8 s -> gg) 27 y mf ( m | + m 2 _ m 2y 



(25) 



x Thc rate (|2"5)l is 8 times smaller than that found by [I] which was used in [UGSGI]. The expression of [7] 
is the correct one, as can be seen by comparing the 1 and 8s cases. For 1 — > gg, diagrams A and A from 
Figure [3] include the color factor 

r^ahfagefehb ~\/% 

where g and h are the adjoint indices of the gluons, while 8s — > gg has the color factor 



^tabc] agej G.hb ^S^Q ac ^bea ^gba fcea) fehb — \/ 5 2 

Squaring and summing over g and h, we obtain Ti^. gg /Ts s ^gg = 4 |?/>i(0)| 2 / |?/>8 S (0)| 2 , as in [7]. The mistake 
of [1] is that they summed over the color index c of the decaying gluinonium in their (2.10). 
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Figure 3: Diagrams for annihilation decays of gluinonia. 



For rrtq ~ trig, these branching ratios are of the order of 5% for a 300 GeV gluino and 0.5% 
for a 1 TeV gluino. There are no other decay modes of the pseudoscalars with considerable 
branching fractions. 

The vector color-octet (8a) gluinonium is not produced by gluon fusion at tree level and 
does not decay to gluons because the sum of diagrams A, A and A' vanishes at threshold. It 
annihilates into qq pairs via diagrams B, B and C in Figure [3] The rate is (based on [7] and 
the m^-dependent prefactor taken from, e.g., the calculation of [32] for qq — > gg) 

r(8 A ^qq) = Y, (4^) ^ l*.(0)| 2 = £ (4^) W ^n s , (26) 

where the summation is over the flavors and chiralities of the squarks. 

Interestingly, whenever a squark has trig = m g its contribution vanishes because of the 
destructive interference of diagrams B and B, which depend on the squark masses rriq, with 
diagram C, which does not. If the squarks are much heavier than the gluinos, only diagram 
C contributes, and (126]) reduces to |7J 

27 

T(8 A -> qq) = —n f a 2 s a 3 s m g , (27) 

where rif is the number of quark flavors. 

Numerically, all the annihilation widths are very small, below 1 GeV. Since these decay 
widths are many orders of magnitude smaller than the gluinonium mass, the gluinonia can 
effectively be treated as 0-width resonances decaying to dijets and ti. The annihilation decay 
modes then have the potential to lead to a precise measurement of the gluino mass. If the 
decay rates are small enough the colored gluinonia will even hadronize with standard model 
particles before annihilating. However, since 1/ao ^> Aqcd, the gluinonium within the hadron 
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will be much smaller than the hadron itself, and so the annihilation decays are insensitive to 
whether or not the gluinonia are confined. 

As expected from (11 ip . in cases when the gluinonium decays primarily by annihilation 
(r ann ^> 2Tg), its width is much smaller than the binding energy. In fact, for the three bound 
states, 

1 27 3 

— r ann (l, 8 S , 8 A ) ~ a 2 s a s (27, — , -n f ) ~ (0.03, 0.003, 0.004) . (28) 
t, h 8 4 

Such gluinonia will appear as distinct states separated from the gluino pair continuum. This 
also provides a self-consistency check on our use of the narrow width approximation instead 
of the full NRQCD formalism which we needed for toponium. 



3.3 Cross sections at the LHC 

As we observed, annihilation decays dominate over other gluinonia decay modes when the 
squarks are all at least as heavy as the gluino. This is of course possible in parts of the 
MSSM, since all the superpartner masses are free parameters. To see whether this parameter 
region corresponds to any viable models motivated by other considerations, we list in Table [1] 
the relevant parameters for a standard set of MSSM benchmark points, the SPS points [10]. In 
many of these scenarios 2Tg 3> r ann so the gluinos will decay much more than annihilate, or, 
if 2Tg > E b , they will even decay before forming bound states, reducing the situation to one 
similar to toponium. These decays will contribute to the usual single-gluino decay signature 
as in the continuum gluino-pair production and the size of this contribution has been studied 
recently in [8]. We are interested in cases where 2Tg < r ann <C Ef, since this is when there is a 
chance to see the annihilation decays. For the pseudoscalar gluinonia (1 and 8s) the rate for 
annihilation decays is significant for the SPS 2 (mSUGRA focus point) and SPS 8 (GMSB with 
neutralino NLSP) scenarios. In these models, two body decays of the gluino are forbidden. 
For the other SPS points, the gluino decays too fast for annihilation decays to be relevant. 
The annihilation rate of the vector gluinonium (8a) turns out to typically have a sizeable 
suppression due to the m^-dependent factor in (|26|) . as also shown in Table [TJ Nevertheless, 
its annihilation is still significant at SPS 2. We can conclude that there are models in which 
annihilation decays may be observable. Since including the effect of a finite branching ratio 
is trivial, in our simulations we will assume Tg = 0. It is this limit that is indicated by the 
last rows in Table [TJ We consider two representative cases: all the squarks degenerate with 
the gluinos (which is roughly the case in most of the SPS points) and all the squarks much 
heavier than the gluinos (like in split SUSY). 

In the narrow width approximation, the production cross sections for gluinonia are simply 
related to their decays through (j57|) : 

<W*) - S (s - Ml) * 2 -^faM 6 (s - M 2 ) , 

^s s (s) - 5 (s - Mi s ) x ^o^l 5 (s - M 2 ) , 
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Table 1: Parameters of various standard MSSM points and the toy scenarios that we will 
analyze. The column (ri™ 1 ) refers to the value without the m^-dependent suppression factor 
in f l26|) . We also compare the gluinonia production cross section cr bound (estimated in the 
narrow- width approximation) with the continuum gluino-pair production cross section <j cont , 
and the last two columns show the annihilation cross sections into dijets and it (while the rest 
of u bound corresponds to single gluino decays like in the continuum). 




a 2 a 3 s 5(s- M 2 ) , (29) 



where the sum in the last line is over squark chiralities. 

Alternatively, the expressions for the cross section can be derived using the near-threshold 
tree-level gluino pair production cross sections, which are (see, e.g., [32] 




Mag -> 99) = 27 ™ 2 S ^o(gigj ->• 99) = t^s_ \- i "^.x ,,u 9 \ 

P 64m? ' (3 6m? ^ L? + ™? I ' 

In the gg channel 1/6 of the production is in the color-singlet representation, 1/3 in the color- 
octet 8s, and 1/2 in 27 (which does not bind), while in the qq channel all the production is 
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in the color-octet 8a [8]. Then, using Eq. (I6"9"j) from Appendix [A], with the wavefunctions in 
Eq. (j!7p computed with the color factors of Eq. (1191) appropriate for each representation, we 
obtain Eq. (j25j) . 

The radial excitations are described by the same expressions as (129]) . but with the annihi- 
lation widths scaled by 

"pann 

pann ^ (31) 

n 6 

In the limit — > the radial excitations can be summed to effectively change r ann by an 
overall prefactor 

^- = C(3)~ 1.202. (32) 

n=l 

For a finite Tg the contribution of the radial excitations is somewhat smaller^] 
The total production cross sections are given by 

a(^^l) = C(3)^3£ ra (M 1 2 )r (33) 
a(pp 8 S ) = C(3)-^£ 55 (M 8 2 g ) r 8s _ gg , (34) 

8s 

a(pp 8 A ) = C(3)|^^£ 9 ,(M 8 2 A )r 8A ^, (35) 



where the luminosities are 



^H^T^t- 1 ' (36) 



where % and j denote the initial state partons and \fs = Ecu is the collider energy. 

These cross sections are shown in Figure H] as a function of the gluinonium mass for the LHC 
at yfi = 14 TeV, 10 TeV and 7 TeV. For this plot, we have used the LO MSTW 2008 parton 
distribution functions [33] evaluated at the scale M = 2m g and neglected both the gluino 
width and the suppression of Ts A ^ qq - due to the finite squark masses. These cross sections 
sum to a few percent of the total gluino pair production cross section, fairly independently of 
the gluino mass over this range. We will discuss this further in Section 16.11 

The cross sections for sample masses rrtg = 300 and 800 GeV, with squarks either degener- 
ate with the gluinos or much heavier, are shown in Table |2j These numbers differ somewhat 
from other results in the literature. Ref. [7] has found cross sections which are a factor of 2 — 3 
smaller than our results. This happens to a large extent because they use instead of 



2 For finite Tg, the excited modes, with their l/n 3 -reduced annihilation rate, still need to compete with the 
same 2Tg. The correct correction can be easily computed numerically. One can also take into account that 
some of the excited states may decay into P states whose annihilation rates are small. This reduces their 
contribution to the annihilation signal. We do not take this possibility into account. Wc also did not take into 
account the direct production of P states and higher angular excitations. 
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Figure 4: Gluinonia production cross section at the LHC as a function of the gluinonium mass 
M ~ 2m g. Solid lines refer to the sum of 1 and 8§ pseudoscalar gluinonia which annihilate into 
gg and it, and dashed lines refer to 8a vector gluinonium which annihilates into qq (including 
it). The three lines in each case correspond, from top to bottom, to LHC center-of-mass 
energy of y/s = 14, 10 and 7 TeV. 



in their expressions for | ■?/>(()) | 2 . The results of [I] for the pseudoscalar gluinonium are about 
4 times bigger than ours (even though they applied certain cuts while we did not). Because 
of the factor-of-8 error in the expression they use for r 8s ^ 99 as we explained in Footnote [TJ 
their combined cross section for 1 and 8s is indeed expected to be by a factor of about 2.4 
larger. However, it also seems that they used a s instead of a s , which should have made their 
result smaller. But they have included a i^-factor of 2.0 which can compensate for that. A 
.fT-factor does need to be included in principle, but its size is still unknown. The effect of 
initial-state radiation on the gluinonium production has been computed in [8], giving a near- 
threshold K-factor that varies from 1.2 for m g = 200 GeV to 1.8 for rrtg = 1 TeV. However, 
this does not yet include the emission of non-collinear gluons which can give an additional 
correction of the same order [8]. Other corrections which can modify the results by ~ 50% 
are higher-order QCD corrections to the potential (pQ), a confining term that may still affect 
the result for the smaller masses, and relativistic corrections suppressed by v 2 . Computations 
of part of the QCD corrections have been recently presented in [9]. Furthermore, we need to 
include corrections also to the annihilation process (final-state radiation etc.). Since some of 
these quantities have not been computed yet, we will not include a K- factor in our simulation 
of the gluinonium signal. We will not include a i^-factor for the QCD background either. 
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Table 2: Gluinonium annihilation cross sections at the LHC. 



4 Gluinonium: simulation 

In the previous section, we found that 300 GeV gluinos can have a cross section for annihilation 
decay as high as 15 pb. Although the cross sections shrink sharply with gluino mass, it is 
certainly possible that the annihilation decays of gluinonia might be seen. The dominant decay 
mode of the gluinonia is to dijets, however the decays to tt may be easier to see because of 
smaller backgrounds. In order to determine what we might learn about the underlying physics 
if these resonances can be found, we now examine the dijet and tt signals and backgrounds 
through Monte Carlo simulation. 

We simulate both the signals and backgrounds at the LHC at 14 TeV using Pythia 
version 8.120 [SUES] with MSTW2008 LO parton distribution functions [33]- For clustering 
the particles into jets, we use the SISCone jet algorithm version 2.0.1 [3S] with radius R — 1 
for the dijet analyzes and R = 0.5 for the tt analysis. 



4.1 Signal and background 

As we saw in the previous section, the bound states of different spin have different production 
and decay modes. The pseudoscalars, which include a color-singlet, $i, and a color-octet, 
$8 S ) couple to both gluons and heavy quarks. In models for which the annihilation decays 
are relevant, the resonances have very small widths. Then we can model them with effective 
interactions of the form 

777 

a = c ) 9 —^ pa ^ G % G % + 4?=r $1 111 q ' (37) 

9 9 

£ 8s = 41 — e^d abc nA G % + 4 s ^ n s Tt > (38) 

frig nig 

where a,b,c and i,j are color indices of the adjoint and fundamental representations, respec- 
tively The factor of m q in the normalization of the qq operators is chosen so that they vanish 
as m q — > 0, as they must by chirality (the rates are approximately proportional to m 2 q : see 
( I2~4j) and (12"5]) ). The vector, V£ , has suppressed couplings to gluons and its coupling to quarks 
can be described by the effective interaction 

^A = 4 A (^)°(n?^- ( 39 ) 
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Figure 5: Dijet search reach (95% confidence level exclusion) at CMS [38] and the gluinonium 

signal with |r7j et | < 1. The steps are due to prescaling. The dashed lines are our extrapolations 
assuming that the prescaling can be reduced by adjusting the triggers. 

All of the coefficients c?; and are fixed by matching to the annihilation rates. One could also 

gg qq j o 

write down operators with additional derivatives, but since in the narrow width approximation 
the momenta of the incoming particles are fixed, we can simply absorb the effect of these 
additional operators into the normalization. 

As a simplifying assumption, since the color matrices do not affect the geometry of the 
processes, at least at leading order, we treat all the gluinonia as color- singlets. This facilitates 
the Monte Carlo simulation. In particular, the pseudoscalar gluinonia reduce to a pseudoscalar 
Higgs with appropriately modified couplings. The vector gluinonium can be simulated as a 
vector Z' with flavor universal couplings of appropriate strength (to the vector current of 
quarks only). We restrict our discussion of the simulation results to the representative cases 



described in Table with m g = 300 GeV (M = 600 GeV) and m g = 800 GeV (M = 



The background for the dijet signals (gg and qq) is QCD dijets. For the it channel back- 
ground, we use the standard model it production, including both gg — >■ it (~ 85%) and 
qq — )■ it (~ 15%) and also single top production (predominantly qq — > tq). We do not include 
fakes from QCD or W/Z + jets since it seems likely the experiments can distinguish these 
backgrounds from it with good efficiency [37], but in any case, top-tagging is a subject largely 
orthogonal to the current work. 
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Figure 6: The invariant mass of the two hardest jets. Left is the 300 GeV gluino and back- 
ground after a pj, > 200 GeV cut, and right is the 800 GeV gluino and background with a 
> 500 GeV cut. 

4.2 Analysis of the dijet channel 

First, let us consider the dijet signatures of annihilation decays. These come from pseu- 
doscalars decaying to gg and vectors decaying to qq. For the m g = 300 GeV case, the jets 
will dominantly have px ^ 200 GeV, and for the m g = 800 GeV case pt > 500 GeV. For the 
light gluino, the 200 GeV pt cut is problematic because it will not survive the high-luminosity 
triggers of either ATLAS or CMS. For example, with an instantaneous luminosity of 10 34 
cm" 2 s -1 , only dijet events with px > 450 GeV can survive prescaling. As pointed out in |38j . 
if we already have information about the gluino mass, for example from cascade decays, than 
the high-level trigger can be adjusted to reduce the prescaling by a factor of up to 1000. 
Therefore, it is still possible that the dijet signatures of the 600 GeV gluinonium annihilation 
decays may be found. In Figure [5j we compare the dijet reach at CMS (for instantaneous 
luminosities of 10 32 , 10 33 and 10 34 cm -2 s _1 ) to the gluinonium cross section. The sawtooth 
shape of the reach is due to prescaling, so we have extrapolated these curves to lower mass, 
assuming the prescaling can be overcome by adjusting the high-level trigger. 

To study the distributions, we simulate the signals and the backgrounds requiring two jets 
with pt > 200 GeV for the m g = 300 GeV case and two jets with p? > 500 GeV for the 
rrig = 800 GeV case. Figure [6] shows the distribution of the invariant mass of the two hardest 
jets in these samples. Note that the signal, whose intrinsic width is a fraction of a GeV, 
broadens to about 5% of the mass of the resonance due to QCD radiation (the scaling of the 
radiation width with the mass of the resonance follows by dimensional analysis, since QCD 
is almost conformal at these energies). The enormous background and trigger issues (for the 
light gluino case) notwithstanding, the gluino mass can, in principle, be determined from the 
location of dijet invariant mass peak. 

Let us now study more closely the behavior of the signal and the background near the 
peak, in the window of 580 GeV < M# < 620 GeV for m g = 300 GeV and 1560 GeV < 
Mjj < 1640 GeV for m g = 800 GeV. In Figure H, we show the distributions of pt and the 
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Figure 7: Distributions of the average transverse momentum of the two hardest jets (top 
row) and the scattering angle in the collision frame (bottom row), in invariant mass range 
580 GeV < Mjj < 620 GeV for the 300 GeV gluino case (left) and in invariant mass range 
1560 GeV < Mjj < 1640 GeV for the 800 GeV gluino (right). 



scattering angle cos 9*. As expected and can be seen from these plots, the angular distribution 
is sensitive to the spin of the resonance. In principle, this information can be used to verify that 
the resonance is indeed the gluinonium. The experiment will measure the sum of contributions 
from two different spins, but as can be seen from Table [2] or Figure HI the spin-0 contribution 
will typically dominate. 

To explore the signal and background rates more quantitatively, consider the 300 GeV 
gluino. In the dijet invariant mass window 580 GeV < Mjj < 620 GeV, including both the 
spin-0 and spin-1 resonance contributions, the signal-to-background ratio is S/B ~ 1 x 10~ 3 
after the cuts pr > 275 GeV, |cos 6 1 *! < 0.5 for both cases from Table El This indicates 
that the signal extraction will be very sensitive to theoretical and experimental systematic 
uncertainties of the background. On the other hand, the statistical uncertainty on the back- 
ground, characterized by S/ \fB, quickly shrinks because we are dealing with a large number 
of events. For example, with ~ 10 fb" 1 there is already a 3a significance. Although the CMS 
analysis [SB] indicates that these signals may be seen (see Figure enhancing the signal-to- 

3 Note however that while the analysis in [38] has taken into account the statistical noise of the background, 
the resolution limitations, and the various systematic uncertainties, it has not simulated the subtraction of the 
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background ratio may be critical to getting a clean sample from which the gluinonium spin 
might be extracted. Possibilities for enhancing the signal might involve jet substructure [3U| . 
color information, additional angular variables, or multivariate techniques. These ideas are 
worth pursuing, but go beyond the scope of the current paper. No matter what, detecting and 
characterizing the dijet annihilation decay signal of gluinonium will be a formidable challenge. 

4.3 Analysis of the tt channel 

Another possible decay mode of gluinonium is to it. As we have shown, this decay rate 
depends on the masses of the squarks as well as the gluino mass. For a 300 GeV gluino and 
300 GeV squarks, the cross section times branching ratio to tt is 780 fb. If the squarks are 
decoupled, the cross section is reduced to 240 fb. Although these cross sections are small, the 
systematic uncertainty on the tt background is much smaller than for dijets, so there is hope 
that with enough integrated luminosity, such resonances may be seen. There has also been a 
fair amount of recent progress in improving the techniques to find boosted tops [101 HU H2, H3] , 
which would come from very heavy gluinonium decays. 

The reconstruction of the invariant mass and the other properties of the decaying resonance 
in the tt channel is more complicated than in the dijet case because the tops decay. As a result, 
the energy and momentum of the resonances are no longer concentrated in two hard jets, but 
in a combination of several jets and possibly leptons and missing energy as well. The cleanest 
channel for tt is the semi-leptonic one, in which one of the tops decays hadronically and the 
other decays to a 6-quark, an electron or muon, and a neutrino. The presence of the lepton 
(unlike in the all-hadronic channel) is very useful for rejecting the QCD background, while 
the fact that there is only a single neutrino (unlike in the all-leptonic channel) allows full 
reconstruction using the on-shell intermediate W. The branching ratio for semi-leptonic tt is 
30%, and, for simplicity, we consider only this channel. The top-tagging algorithm we use is 
described in Appendix O 

The resulting it invariant mass peak is fairly sharp, as shown in Figure |HJ To study the 
signal in more detail we now restrict to the mass window 575 GeV < M t j < 625 GeV. The 
Pt and cos 9* distributions in this window are shown in Figure [9j The difference between 
the angular distributions of the scalar and the vector is much smaller here than in the dijet 
case because the top is not much lighter than the 300 GeV gluino. In the limit of very heavy 
gluinonium, the differences between the spins reappear, but the cross section for the tt signal 
drops too fast to make this channel useful (from Table [2, for m- g = 800 GeV, cr t j ~ 3 fb if 
the squarks are decoupled and an order of magnitude smaller if they are degenerate with the 
gluinos). 

For the case rrig = rrig = 300 GeV we obtain S/B = 1.6 x 10 -2 and S/yfB = 3 at 500 
fb -1 after applying the cuts px > 200 GeV, |cos#*| < 0.5. These estimates are somewhat 
less promising than what has been found in generic analysis of resonances decaying into tt by 
ATLAS [44] . For example, they find that 300 fb" 1 of luminosity will allow a 5cr discovery of 
a 600 GeV resonance if its cross section is above ~ 700 fb, which is approximately the cross 
section we have. In contrast to the dijet channel, the it search is entirely statistics limited. 

background from the data, which we guess is likely to be a critical issue with the small S/B that we have. 
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Figure 9: Distributions of the average pt (left) and the scattering angle of the it pair in the 
collision frame (right) for the 300 GeV gluino in invariant mass range 575 GeV < M t j < 
625 GeV. 



The reach could be improved somewhat by inclusion of other it decay channels, and through 
an optimization of the cuts and top-tagging algorithm. We conclude that it is a viable, but 
high luminosity, channel in which gluinonium, if it exists, could almost certainly eventually 
be found. 

4.4 Estimate for the diphoton channel 

Since the gluinos are not charged under the electroweak group, the gluinonia can decay into 
electroweak bosons such as 77 only through loops, with branching fractions of ~ 10 -5 [9]. 
Nevertheless, it might still be useful to consider the diphoton channel because of its much 
smaller background compared to the dijet and it channels. We have not performed a complete 
analysis, but it is easy to get a rough idea by comparing to the analysis of the 77 signal of the 
squarkonium in |35] which we will mention in the next section. It turns out that the small 77 
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branching ratio for the gluinonium compared to the squarkonium is roughly compensated by 
the larger color factors in the gluinonium production process, and the signal may be observable 
if the gluino is relatively light Jj For rrig ~ 300 GeV, we will have S/B of about 10%, which is 
much better than in the dijet or it channel, but 3a significance will require almost 10 3 fhr 1 of 
luminosity 

5 Squarkonium and other QCD bound states 

So far we discussed the bound states of color-fundamental fermions (toponium) and color- 
adjoint fermions (gluinonium). We have seen that relative to the decay rate of the tops, 
the annihilation decays of toponium have a too small branching ratio to be detectable. On 
the other hand, the gluino is sufficiently long-lived in certain regions of parameter space, 
so gluinonium may have observable decays to dijets and it. Let us now discuss the other 
possible bound states of particles charged under QCD. In the MSSM, the only other colored 
particles are squarks, which are scalars in the fundamental representation. There exist strong 
production processes (see for example [I2JE2]) for a squark-antisquark pair 

99 -> QiQi , QiQj -> Mi , (40) 

a pair of squarks 

qiqj -> Mj , (41) 

or a squark and a gluino 

Qi9 -> Qi9 • (42) 
Let us discuss the possible bound states of these particles. 

5.1 Squarkonium 

Squarkonium is a bound state of a squark and an antisquark. As in the case of toponium, the 
color decomposes as 

3 g> 3 = 1 © 8 , (43) 
and only the color singlet can bind. The color factor for the singlet bound state is 

C x = | . (44) 

In most scenarios squarkonia do not form because the squarks decay too fast. For the 
gluino, we could make the single-particle decay rate arbitrarily small by making all the squarks 
heavy. However, because squarks have electroweak quantum numbers and are heavier than 
the lightest supersymmetric particle (LSP), they can generically decay as 

q — > q' W or q — > qx , (45) 
4 We thank Cliff Burgess for suggesting us to check this possibility. 
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Figure 10: Diagrams for the dominant annihilation decay of stoponium. The first two diagrams 
do not actually contribute because they vanish at threshold. 



Model 


m~ t (GeV) 


E b (GeV) 


2r,-(GeV) 


T ann (GeV) 


SPS 5 (mSUGRA with light stop) 
All other SPS points 


262 
> 400 


2.3 

~ 5 


0.087 

> 3 


0.0028 
~ 0.004 



Table 3: The mass of the lighter stop m t -, the binding energy Eb, twice the stop decay rate 
2r t ~, and the stoponium annihilation rate r ann , for the SPS benchmark points [TU]. In all the 
SPS scenarios, the stop decays roughly equally via t — > x°t and t —> y + 6 (only the latter is 
present in SPS 5). 



where x is a chargino or neutralino. The rate for these decays is in general comparable to the 
binding energy 

Eb = ^ (46) 

of the would-be squarkonium and we obtain a situation similar to the toponium in which the 
bound state is very broad and it decays primarily by single-particle decays rather than the 
annihilation decays. 

The only possible exception is stoponium. There, the stop can be lighter than the top 
plus LSP, so the weak decays are forbidden. This is not the case in most of the parameter 
space of the mSUGRA ansatz. In SPS 5, even though the width is smaller than the binding 
energy, single-stop decays still dominate over the annihilation decays (see Table [3]). However, 
the situation is much more favorable in certain other motivated scenarios discussed in more 
detail in [451147]. 

Stoponium (J PC = ++ ) is produced primarily by gluon fusion. Its annihilation into gg 
via the diagrams in Figure [10] has the rate 

r(l ^99) = ^§ \m I 2 = §«^m 9 -. (47) 

Unfortunately, this is about 200 times smaller than for gluinonia of the same mass mostly due 
to the different color factors (compare to fT22l) . and note that in the gluinonium case there are 
also contributions from the octets), and twice smaller than for a heavy quarkonium of the same 
mass (see (jlOp). and the production cross section, which is given by (1331 . is correspondingly 
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smaller. As a result, the stoponium cannot be seen in the dijet or tt channels even when the 
annihilation decays dominate. 

However, since the squarks are charged, the annihilation of stoponium to 77 is relevant. 
Although the branching ratio to 77 is only ~ (8/9)(a;/a s ) 2 ~ 0.005, for certain regions of 
MSSM parameter space the cross section may still be of order 1 fb which makes it possibly 
observable. Other stoponium decay modes including ZZ and W + W~ could possibly be visible 
as well. Stoponium decays have been studied in detail in [3S1 HH1 SSI EDI EI] (note that [1H1 EH] 
contain certain factor-of-2 errors - see [5~2"l ESI US])- 

5.2 Di-squark and squark-gluino bound states 

Instead of a squark and anti-squark, another possible bound state might be two squarks or a 
squark and a gluino. For di-squarks, the relevant color decomposition is 

3 <g> 3 = 3 © 6 (48) 

and the color factors are 

so the triplet has an attractive potential and can potentially form a bound state. For a 
squark-gluino pair, the decomposition is 

3©8 = 3©6©15 (50) 

and 

3 1 1 

C3 = 2 , Cq = - , C15 = -- , (51) 

which means that the potential for the triplet and the sextet is attractive and could potentially 
support spin-| bound states of these two particles. 

However, the di-squarks and the squark-gluino bound states suffer from the same problem 
as the squarkonium, namely the too rapid decays of the squarks. Furthermore, even the special 
cases in which the lighter stop is sufficiently stable are not available here because, unlike in 
the case of the stoponium, producing a di-stop or a stop-gluino bound state via f l4"Tj) or f l4"2"j) . 
respectively, would require having top quarks as incoming partons. 

5.3 Other QCD bound states 

There are many more possible bound states of hypothetical new colored particles. For example, 
superheavy quarkonia have been studied in [5^1 [55] and bound states of 5'L r (2) Z/ -doublet color- 
octet scalars in [56]. It would also be interesting look at bound states of Kaluza-Klein (KK) 
excitations of quarks and gluons in theories with extra dimensions. For example, Universal 
Extra Dimension (UED) models have TeV-scale flat extra dimensions with all the standard 
model fields propagating in the bulk, so they all have massive excitations. For certain regions 
of parameter space KK quarks can be sufficiently stable to form bound states [571 158] - Even 
though the KK particles are color adjoints or fundamentals, like the gluino or squarks we have 
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already considered, the details of the phenomenology will be somewhat different because of the 
different spins as well as the different range of possibilities for the single-particle lifetimes, and 
it would be interesting to study this in more detail in a future work. If annihilation decays can 
be detected at the LHC, one interesting application would be to try to use them to distinguish 
extra dimensions from supersymmetry. From Figure 0, there is clearly some spin information 
contained in the dijets from annihilation, but isolating a clean sample to see this effect may 
be impossible in practice. The size of the cross sections for the different annihilation channels 
might also be a useful discriminator, when combined with other information. 

6 Bound states with new forces 

The previous sections have all looked at bound states for which the binding force is the strong 
force of QCD. We have seen that being charged under the strong force allows particles to be 
produced with substantial cross sections, and substantial binding and annihilation rates, as 
high as tens of picobarns. Our discussion of the phenomenology of colored bound states was 
centered around supersymmetry, in particular, the MSSM, which has new particles but no 
new forces as compared to the standard model. However, there are plenty of models which do 
have new forces. Technicolor is an obvious example, but even within supersymmetric models, 
dynamical supersymmetry breaking often involves something beyond 577(3) x SU(2) x U(l). 
It is certainly possible for there to be a new force, even stronger than QCD, which would lead 
to quite substantial annihilation decay rates. 

Let us suppose there is a new force associated with a new SU(N) gauge group. If the 
associated strong coupling scale A^r is sufficiently large (say 10 TeV), then all the particles 
produced at the LHC would be confined into S77(AQ-singlets. This case is not interesting 
for the present paper, as there is no way to distinguish particles bound by this force from 
fundamental resonances. Instead, we are interested in a situation in which the LHC will be 
able to produce particles with mass m ^ A N > Aq CD . We further assume they are charged 
under both QCD and SU(N), so that they will also have substantial production rates. Recall 
that the annihilation rates are proportional to |?/>(0)| 2 oc a^. Since these rates scale as the 
third power of the coupling, the new gauge coupling does not have to be much larger than the 
QCD coupling to drastically increase the rate. 

To distinguish annihilation decays from simple resonances, we also need to have some other 
way to observe these new particles. For example, this may be through single particle decays, 
as with the gluino. In fact, there are many models in which weird things happen when there 
are new forces, such as displaced vertices in hidden valley models [59], or macroscopic strings 
in quirks [60]. We will focus only on annihilation decays, leaving their coordination with the 
rich field of new-force model signatures for future consideration. 

6.1 Increasing the annihilation rate 

For gluinos of phenomenologically relevant masses, we saw that as the gluino becomes more 
stable the gluinonium annihilation cross section saturates at a few percent of the continuum 
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gluino pair production cross section. Let us discuss which factors set the scale for this number 
and what can make it larger, and then consider the situation in which the binding particles 
have an additional attractive force, beyond the SU(3) of QCD. 
The bound state production cross section is given by 

jb ound ^ £(4m 2 ) /".bounds j./;. C{Am 2 ) a 2 [JJn ^ 2 r( A _^a 2 C 3 a 3 s 



m 



[ <x bound (s) dVS ~ ^ |^(0)| 2 ~ £(4m 2 )^^ , (52) 

J m m 4 m z 

where C(M 2 ) is the parton luminosity, Eq. ( )36l) . If we neglect the dependence of the parton 
luminosities on the energy, then the parton-level continuum production cross section <T cont (s) ~ 
a 2 /s gives 

a cont „ ^m 2 ) r ^ con ^ ^ ^ «| 
ml m 2 



so that 

abound 



- ~ C d c^ . (54) 

v 

There are additional factors coming from the fact that only part of the color configurations 
bind, and the fall off of the parton luminosities with energy, but this qualitative scaling still 
holds. The point is that because C and a s are raised to the third power in this relation, 
the fraction of annihilation decays can be significantly enhanced over the continuum if the 
coupling or color factors are larger than they are for QCD. 

It is interesting to also consider how large the fraction of annihilation decays can be, while 
the binding force is still perturbative. When the coupling a s is increased, it not only increases 
the binding rate, but also adds a radiative correction to the continuum production. This can 
lead to a significant Sommerfeld enhancement of the total cross section. From Eq. (I72p . the 
additive correction to the continuum production near the threshold is roughly 

a 2 

Aa cont (s)~^-Ca s . (55) 
m z 

This extra contribution needs to be integrated up to roughly the binding energy E^ ~ C 2 a 2 m. 
This gives 

A^cont i /•2m+E b bound 

=?— t — / Aa cont (s) dV§ ~ C 3 a 3 s ~ ^— . (56) 

^,-corit ^cont / v ' A ^cont v ' 

Thus for large Ca s the continuum rate will also be proportionally larger, and there will never 
be more than ~ 50% annihilation decays (with a s < 1). 

In the more extreme limit in which the theory becomes strongly coupled and confining, 
the bound states will look just like simple resonances decaying either back into the standard 
model or within the strongly-coupled sector itself. But, in the regime of interest, when a s is 
perturbative but (Ca s ) 3 is larger than for gluinonium, the binding particles will be observable 
on their own and the annihilation decays of their bound states will be significantly enhanced. 
In the next subsection we exemplify this with a toy model in which the particles are charged 
under an additional force. 
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6.2 Example model 

Suppose there exist color adjoint fermions, like the gluino, which also transform under some 
representation R of a new SU(N) with a coupling constant a at. Their pair production pro- 
cesses will be like for the gluino, but without the contribution from virtual squarks. For qq 
initial states, only the anti-symmetric color adjoint state 8a is produced near the threshold. 
For gg initial states, pairs in the 1, 8s and 27 color configurations will be produced with 
relative cross sections o" (l) : <t (8s) : ^o(27) = 1:2:3 [S] • Since the initial states are SU(N) 
singlets, only SU( iV)-singlet configurations can be produced. Because of the sum over final 
states the rate is enhanced by an additional factor of the dimension of the representation, Dr. 

For bound states, the cross-sections will also include the appropriate |?/;(0)| 2 factor for each 
color configuration. For computing those, note that the couplings in the potentials (CQ) are 
given by the simple replacements 

Ca s 4Ca s + C^ajy , (57) 

where C is the relevant QCD color factor from (fl9|) for the 1, 8s, 8a or 27 states, and Cr is 
the quadratic Casimir of the representation R of SU(N). For example, for the fundamental 
representation of SU(N), Dr = N and Cr = (N 2 — 1)/2N, and for the adjoint representation 
D R = N 2 -1 and C R = N. 

To be explicit, the cross sections for the pseudoscalar bound states in the (1, 1), (8s, 1) and 
(27, 1) representations of SU(3)qct> X SU (N) will be like the cross section of the 1 gluinonium 
of the MSSM times the factors 

D / 3*+fW |' 2Dr / 3/2B . + C R a N V f -a, + C R a N V 

\ 3a s J \ 3a s J \ 3a s J 

respectively. The (8a, 1) vector will be like the 8a gluinonium times 



3/2 a s + C R a N 



H 3/ 25 ; ) ■ < 59 > 

Notice that even in relatively small representations and with order 1 values of a^/Hs the 
annihilation decay rates will be orders of magnitude larger than anything in the MSSM. 

Figure [11] compares the dijet rates through annihilation decay to the CMS dijet search 
reach for several toy models, including the one with the MSSM gluino. Even when the addi- 
tional force has the same strength as QCD (a at = a s ), the cross section is enhanced by more 
than a factor of 10 if our "gluinos" are charged under the fundamental representation of a 
new SU(3) group. If they are charged under the adjoint of a new S77(4), the enhancement 
is almost 1000. If the new force is somewhat stronger, a N = 3a s (which corresponds to a 
confinement scale of about ~ 10 GeV, assuming for simplicity the same running as in 
QCD), the signal will be 10000 times larger than the annihilation decay signal from MSSM 
gluinonia. 

The bound states may also annihilate into hidden sector dijets. But note that if the QCD 
and hidden sector couplings are comparable, the branching ratio into QCD dijets will still be 
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Figure 11: Dijet search reach (95% CL exclusion and 5a discovery) with 100 pb _1 at CMS [55] . 
the MSSM gluinonium signal, and several models in which the particles are bound also by 
new forces. 

of order 1, and the dijet rate will still have an enormous enhancement over annihilation decay 
rates in the MSSM. The phenomenology of decays into the hidden sector is very interesting, but 
model dependent, so we will not consider it here. In summary, there is plenty of parameter 
space for complete models of this kind which are consistent with the various experimental 
bounds, in which such annihilation decays can provide spectacular signals at the LHC. 

7 Summary and conclusions 

In this paper we have considered the annihilation decays of a variety of possible bound state 
particles at the LHC. A number of examples were discussed at length, including bound states 
of top quarks, toponium, and bound states of gluinos, gluinonium. We also discussed bound 
states of squarks, squarkonium, and bound states associated with possible new strong forces. 
We have discussed how the spins and color representations of the particles determine the 
possible bound states, what determines their production channels and their cross sections, and 
how the mass spectrum of the model determines whether the bound states form and whether 
they decay by annihilation or through single-particle decays. Whenever the observation of 
the annihilation signal (or its absence) is feasible, it can provide valuable direct or indirect 
information about the underlying particles. 

In the first and only standard model case, toponium, we found that the width of the 
top quark itself precludes any reasonable hope of finding toponium. Not only does the top's 



27 



relatively large width force the top quark to decay much faster than it can annihilate, but it 
even decays too fast for a well-defined bound state to form, broadening the annihilation decay 
spectrum. There is no sharp distinction between the dijets coming from toponium annihilation 
decays and the contribution that virtual top quarks give to continuum dijet production. One 
might have also hoped to find evidence for toponium from the detailed spectrum of the weak 
decay products, but since there is no sharp feature at the it spectrum, the resolution required 
is beyond the ability of the LHC. 

Although the top lifetime is fixed, gluinos may live long enough for the gluinonium to 
decay predominantly by annihilation in certain regions of parameter space, including some 
of the standard SPS points. We found that gluinonia may have substantial decay rates both 
to dijets and to it, and that both of these decay modes, although difficult to see, may be 
detectable at the LHC if the gluino mass is of the order of several hundred GeV. Looking at the 
mass distributions we found that both channels could provide peaks for mass determination. 
Moreover, the angular distributions in the dijet channel may help determine the spins of the 
resonances and thereby put constraints on the properties of the binding particle. However, 
due to the low signal-to-background ratio it will be challenging to extract spin information 
from the data. 

For the case of bound squarks, squarkonium, observing annihilation decays into either dijets 
or it pairs is practically impossible. First of all, despite the fact that the production is through 
gluon fusion, squarkonia cross sections are always several orders of magnitude smaller than 
gluinonia cross sections because of color factors. Moreover, for generic models, the squarks, 
like the top, will decay too fast relative to the rate of the annihilation processes or even before 
they form bound states, and the large width will smear the bound state resonance signal. Only 
stops may have a long enough lifetime so that the annihilation decay signal might possibly 
be seen. For these the best hope is using the small branching ratio to annihilate into two 
photons. Although the cross section for these decays is less than 1 fb for stop masses above 
200 GeV, this channel may be feasible for high luminosity. In fact, this channel might even 
be feasible for the gluinonium because of its large color factors even though the gluinos can 
annihilate into two photons only through loop diagrams. We also considered bound states of 
a squark and a gluino and two squarks (as opposed to a squark-anti-squark). Neither of these 
have visible annihilation decay signals. 

Finally, we looked at annihilation decay cross sections in models with new forces. With a 
new confining force, the rate for annihilation decays may easily be much larger than in the 
supersymmetric models which only bind with QCD. There is a huge space of possible models, 
and just restricting to one simple class, we found that the signals for decays of these bound 
states are easy to see. Of course, in the limit that the new force is infinitely strong, the 
annihilation decays are just decays of a new resonance. However, in intermediate regimes it 
is easy to have both annihilation decays and other signatures, so that the annihilation decays 
provide complementary information. 

Overall, it seems that annihilation decays of particles bound through QCD are difficult, 
but possible, to see, while particles bound through stronger forces will show up quickly. If 
the binding is not too strong, as with QCD, we will almost certainly have information about 
other decays of the constituent particles before we find the bound states. This will be useful 
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in directing our search - knowing the mass ahead of time can improve the search strategy, 
especially in the dijet channel. Conversely, if the binding is very strong, exploring what the 
resonance may be a bound state of may help direct searches for other physics associated with 
the binding particles. In general, annihilation decays of bound states have the potential to 
provide an additional handle on the new physics we are likely to see at the LHC. 
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A Bound state formalism 



For computing the bound-state effects, as well as finite-width effects, it is useful to consider 
the Green's function G(x, E) of the Schrodinger equation describing the two-particle systemQ 
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+ V{r) 
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G(x,E) = 5 (3) (x) 



(60) 



where E = \/§ — 2m. For zero angular momentum, the solution is 
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where F{x) = T'(x)/T(x) is the digamma function, 7^ = — ^(1) ~ 0.577 is Euler's constant, 
C is the color factor (|2]) and E n = C 2 a 2 s m/An 2 are the energies of the radial excitations. The 
solution describes both the bound states and the near-threshold continuum. Particles with 
width T (which gives width 2r to the bound states) are described by E — > E + iT p3J [6TJ . 
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A.l Production cross section 

By the optical theorem, the production cross section will be proportional to ImG(0, E), thus 
it can be written as (E3 ESI E2 El EQ 

A /rn , lmG(0,E) . . 



5 Here x is the distance between the particles and we will be interested in G(x = 0, E) because we want 
to create the particles at zero separation and annihilate them at zero separation at a later time (the time is 
represented by its Fourier variable E). 
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where <5"o is production cross section of the pair of particles (in the particular angular mo- 
mentum and color representation) without the binding corrections and Go is the free Green's 
function, that is flBTj) with C = (and T = 0), namely 



m^ / K / F] 

G (0,E) = --J , ImG (0,£) = —J -6(E) ~ (63) 

An V m An V m An 



where /3 = y 1 — Am? /s is the velocity of the particles in their center-of-mass frame. Then 

where we wrote the last expression in terms of the 2-particle phase space 

*2= [ , ^ o f , d > o M V (p 2 - Pl ) 5 (mf ^ A _ A\ = L (65 ) 
J (2nf2p\J (2vr) 3 2p° v ; VH P ; V 2m 2m J 8n v ; 

since the combination <9"o/$2 remains well-behaved as we go below threshold. 

For an attractive potential (C > 0), the last term of f lBTj) grows large near the bound state 
energies E = —E n if E n 3> T: 

G(0E)~-Y ^ ( ° )|2 ~-V 2M "l^°)l 2 (66 ) 

p AM 2 Y 

ImG(0, B) = £ 1^,(0)1' (g + gi)3 + p - £ l^(0)P (| _ M32 " +4My2 , (67) 

where 1^(0) | 2 = (7 3 a 3 m 3 /(87m 3 ) refers to the hydrogen-like wavefunctions describing the 
bound states and M n = 2m — E n are their masses. In the limit r — > 0, 

Im G{0, E) = irJ2 \M0)\ 2 8{E + E n ) = ^ M n |^ n (0) | 2 $x , (68) 

n n 

where $1 = 27r5(s — M 2 ) is the usual single-particle phase space. The last expression makes 
the agreement of the whole formalism with f fl5|) manifest^) However, (fT5j) is more general since 
it holds also for non-Coulombic potentials (but it is valid only for narrow bound states). 
More generally, the imaginary part of (ISTj) is [621 [HI [6] 
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ImG(0,£) 



^_ _x /«+ \ A „ Ca s T/n + 2i;+ (VF+P + E n ) 

v + + Ca s tan 1 — + > E n - '- 

W ti (E + E n ) 2 + r* 



(70) 



6 In the case that the two particles are identical, $2 in (|65[) will actually be twice the phase space. This is 
consistent with the factor of 1/2 discussed after (TT6|) . 
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where 



" ± = ^ 2^ ' (71) 



To see what (1701) means, suppose r — > 0. For _E < 0, Eq. ( 1701) reduces to a sum of narrow 
Breit-Wigners ( 1671) in the attractive case and vanishes in the repulsive case. For E > 0: 

ImG(0,£) = ^^- 1 ^_ (72) 

4 1 - exp (-TrCa s /l3) 



For a repulsive potential, the imaginary part (jT2j) is suppressed to zero at the naive threshold 
E — 0, but for an attractive potential it starts with a finite value 

ImG(0,£^0 + ) = ^^. (73) 

at the threshold. 

A. 2 Annihilation processes 

Consider the process 

ab ^ a(3 ^ AB . (74) 
When a(3 form a narrow bound state, the standard resonance production formula (1951) gives 

(7 a6 ^ M) ^B(s) = — — 2 [ X2 if a = b } , (75) 

^>a^b [s — M l ) + sT z 

where M ~ 2m is the mass and J is the spin of the bound state and D x denotes the dimension 
of the color representation of the particle x. 

When the bound state is not narrow or we want to obtain the contribution from the near- 
threshold a/3 continuum as well, we can use the fact that G(0, E) is the two-point function 
describing the pair, so the matrix element (squared) is proportional to \G(0, E)\ 2 , while all 
the short- distance factors can be taken from the narrow bound state case. Dividing (1751) by 
the narrow-width \G(0,E)\ 2 from (I6"6jr 1 and multiplying by the general ^(O,^)) 2 we obtain 

a ab ^ AB (E) _ 4m 2|^ ( o)|4 \°^ E ^\ [ X2 if o - 6] . (76) 

Note that |^(0)| 4 in the denominator cancels the |^(0)| 2 factors in the annihilation rates, so 
the prefactor of \G(0,E)\ 2 in (JTB"]) depends only on the short-distance physics. We can see 



7 More specifically, \G(0,E)\ 2 from (|66p includes the interference with the radial excitations, which we 
should have, in principle, included in (|75[) as well. We ignore this for simplicity of presentation. 
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this more explicitly by comparing (169 p and (l9~Tjn which allows us to express the bound state 
annihilation rates in terms of production cross sections of free particles, giving 



Cab^aB^AB 
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(77) 



Written in this form, the expression is useful also for color representations that do not bind. 
For the bound states in the narrow- width limit, f|66|) and f )67|) show that 
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where in the last step we promoted the decay rate to include not just the constituent particle 
widths T but also the (sum of all the) annihilation rates r ann : 
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This makes sense because then (ITS"]) or ([77]) reduces to a branching ratio times the production 
cross section (|6T 
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When the constituent particle width T is not much smaller than the binding energy Eb, 
\G(0,E)\ 2 is no longer proportional to ImG(0, E) so the branching ratio is not simply de- 
termined by r^^^/Ttot- In fact, the annihilation rates are not even well-defined in this 
situation because the bound states are not well-defined. Nevertheless, we can still use (1771) . 

However, this method cannot be used when the annihilation rate r ann is the one that 
becomes comparable to Ef,, because the width T in this formalism is the single-particle decay 
width, without taking into account the possibility of annihilation (see, e.g., the derivations 
in [HI [61]). But this regime will not usually be relevant to actual systems because of ( |TTi) . 

Let us now study the behavior in the continuum. Note that for \E\,T < E b we can 
approximate (16"T]) by 
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3 This expression applies if A and B are massless particles with 2 polarizations, but can be easily generalized. 
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Using (I77p with (I83p and (1M|) with (173"]) . we obtain the branching ratio of the annihilation-like 
processes in the near-threshold continuum: 
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With typical strong- interaction diagrams we will have a^ B ^ a ^ /$2 ~ a 2 /m 2 giving 
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A. 3 The case of tt 



Let us apply the results of this Appendix to computing the curves in Figure [U The near- 
threshold tree-level tt production cross sections in the various spin and color channels are 
(see, e.g., [65]) 

\a^ ls ^ , a^ 1SoW , 3Sl(8) ) = (1 , J- , I\ x ^ ■ (86) 
I ' ' J \ 96' 192' 9 J K ' 

From (El]), the production cross sections for these three channels are 

4vr 2 o- 2 f 1 5 1 1 

^ = -^ix|-ImG 1 (O ) £0 ) — ImG 8 (0,£), - ImG 8 (0, £) j , (87) 

and from (1771) the annihilation cross sections (for A,B = a, b) are 

- ann (£) = ^ >< ^(O,^ 2 , H |G 8 (0,E)| 2 , | |G 8 (0,£)| 2 } . (88) 

In practice, when the potential is repulsive \G(0, E) \ 2 is highly suppressed in the near-threshold 
region, so only the color-singlet contributes to the annihilation signal in the case of tt pairs. 
For the near-threshold continuum, just for the color- singlet, ( 184"]) gives the branching ratio 

A 2 2 

a ^ a ^ ~ 7g + 71 Cr^a 2 . = 2.4 x 10" 3 . (89) 

B Resonance cross sections from decay rates 

Consider the process 

a b -> M r A £ . . . , (90) 

in which partons a and 6 (quarks or gluons) produce a resonance of mass M, width T, spin 
J, in a specific color representation, that then decays into particles A and B (or more). The 
resonance propagator contributes the denominator 

1 
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The tensor structure of the propagator can be written as a sum of products of polarization 
vectors which will be absorbed in the decay rates that we will now discuss. The final part of 
the process (1901) is the same as in the decay M — > AB (by rotational invariance, the total 
decay rate is the same for every polarization of M) , so it contributes 

2MF M ^ AB . (92) 

Similarly, the initial part of the squared matrix element, summed (for the purpose of the 
derivation) over all possible directions for the incoming particles and over the polarizations of 
a, b and M, is proportional to T^-mb- It is given by 

(2J+1)^ 4tt 

2D a -2D b s [ 1 y ' 

where we undid the integration over the different directions (which is trivial after summing 
over the polarizations), denoted the dimension of the color representation of X by Dx, and 
assumed a and b to have two helicities. This gives 

Air (2 J + 1) D M M 2 T M _> ab T M ^ AB 

a ab ^M^AB\s) = — — o x2 it a = b . 94 

1 ' D a D b s (3 - M 2 ) + M 2 T 2 [ 

We should note, however, that the actual mass of the resonance is \/I rather than M, and we 
need to make this replacement throughout the calculation: 
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D a D b (s-M 2 ) 2 + §T 2 



0-ab^M~>AB{S) = ^ — 2 ^ n [ x2 if a = b] , (95) 



where the widths TM->ab, Tm->ab and T should be evaluated assuming resonance mass Vs 
instead of M. In the narrow-width limit T <C M, fl95|) becomes 

/~\ 27T (2J + 1) D M YM^ab^M^AB n , /r 2\ r o -c L l fnr\ 

a ab ^ M ^AB{s) ^ =-= — 2ttS(s - M ) [ x2 if a = b ] , (96) 

D a D b M 1 

and if we further sum over all the possible final states in (|9~U|) this simplifies to 

/^\ 2-7T (2 J + 1) D M r^f^ob 2 r 

o-ab^M(s) ^ =r-=: -— 27r5(s-M) [ x2 if a = 6]. (97) 

U a L) b M 

As expected, the production cross section is independent of the available decay modes. 



C Algorithm for tt reconstruction 

This is the algorithm we use for the reconstruction of the tt pair in the semileptonic channel 
with an electron or a muon. It includes the most essential features of the algorithms studied by 
ATLAS [50] and CMS [37]. The efficiencies given below were obtained for signals of scalar 
and vector 600 GeV resonances and the standard model background that was described in 
Section I4TT1 The efficiencies are listed in the order scalar/vector/background. 
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1. We select events with exactly one electron or muon with pt > 35 GeV and rapidity 
\y\ < 2.5. Efficiency: 30%/28%/22%. 

2. We save all jets with p T > 20 GeV, but require to have at least 4 jets with p T > 30 GeV. 
Efficiency: 75%/57%/64%. 

3. We require that the 4 hardest jets have rapidities \y\ < 2.5. Efficiency: 75%/75%/67%. 

4. We start with reconstructing the leptonically decaying top. We compute the transverse 
momentum of the neutrino from the transverse momenta of the lepton and all the jets. 
By requiring W to be on-shell, we obtain a quadratic equation for the longitudinal 
momentum of the neutrino, and keep both solutions. If the discriminant is negative, we 
try reducing neutrino's transverse momentum by up to 15% till the solution is possible. 
Efficiency: 76%/79%/76% (neutrino momentum reduction was used in 7%/7%/8% of 
these cases). 

5. We then find which of the jets best reproduces the mass of the top when combined with 
the W. We choose the neutrino momentum assignment that gives the better top mass 
accuracy. We require top mass accuracy of at least 10%. Efficiency: 87%/87%/88%. 

6. To reconstruct the hadronically decaying top, we take the remaining jets and find the 
assignment of three of them that reconstructs the W and the top most successfully, 
based on the criterion of minimizing the quantity 



Am w \ 2 / Am t \ 2 
0.1m w J + \0.2m t ) 
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where Amw and Am t are the differences between the reconstructed and the actual 
masses. We require that the best assig nment has e = e bes t < (efficiency: 47%/41%/46%) 
and the next-to-best assignment has e > 3eb es t (efficiency: 24%/34%/26%). 

7. To eliminate many of the events that were not reconstructed correctly, we require that 
the azimuthal angle between the two tops is at least 160° (efficiency: 48%/67%/46%). 

The fraction of events that remains at the end of the algorithm is 0.6%/0.8%/0.34% out of all 
the it events. 
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